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Abstract 



In this work, we calculate the solutions of the Rarita-Schwinger equation with 
^ ■ the inclusion of the electromagnetic interaction. Our gauge and coupling pre- 

scription choices lead to Dirac-type solutions. One of the consequences of our 
results is that all inconsistencies related to noncovariance and noncausality 
are avoided and the Landau level occupation of particles turns up to be quite 
different from the usual spin 1/2 particle system occupation numbers. 
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1. Introduction 

It is well known that in heavy ion collisions, when dense matter is formed, 
it is rich in delta resonances Even if low energy physics is considered, 
this kind of resonance is very important. In the irN interactions for example, 
a problem that has been studied for a long time, the A resonance plays a 
central role 0, S B S|- The scattering amplitude, saturated by the A, is 
dominated by this contribution in the n + p channel and it is also important 
in the other ones. The importance of spin 3/2 baryons is also observed for 
other hadron interactions, as in hyperon interactions, where resonances with 
spin 3/2 (S*(1385), S*(1533) and others) [6] also dominate the scattering 
amplitudes in many processes. Studies involving the chiral scalar form factors 
of the nucleons [5| and hyperons also show the dominance of the A and 
hyperon resonances. 
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Hadronic interactions are also important in high energy physics. When 
considering particles produced in a close region, they may be considered as 
comovers in an expanding system, and their relative energy is small. There- 
fore, the effect of the resonances may have a relevant contribution in the 
whole process. 

The above mentioned facts were a strong motivation for the delta baryons 
to be considered as possible constituents of neutron stars 

BUB EH. The 

appearance of particles in neutron stars described by relativistic models de- 
pends on the symmetry energy density dependence [jjj], the choice of cou- 
plings between baryons and interacting mesons (still unknown) 0,0] and the 
possible appearance of pion and kaon condensations 0, 13, HI]. In [i| the con- 
tribution of the four delta particles (A~, A , A + ,A ++ ) is taken into account 
in a mean field theory (MFT) approach and the related Rarita-Schwinger 
(RS) equation of motion 0] is solved self-consistently with the other equa- 
tions of motion that describe the other baryons and mesons. The population 
of A particles in neutron stars for different asymmetries (/9-equilibrium is not 
enforced) is then obtained. In 1(|, the influence of A resonances on the max- 
imum star mass is investigated with two different relativistic models. Indeed, 
the inclusion of A baryons leads to an increase of the maximum mass within 
the non-linear Walecka model [15j and a decrease within a chiral hadronic 



model 16 



Neutron stars with As are revisited in ll[ and they seem to 
suppress the hyperon abundances considerably. In the last years, neutron 



17 



stars subject to strong magnetic fields, named magnetars by Duncan 
have been a source of intense investigation 18|. If one assumes that A par 



tides have to be seriously treated as possible constituents of these compact 
objects, the RS equation with the inclusion of the magnetic field 
has to be taken into account. 

Recently, the influence of magnetic fields on the QCD phase diagram 
for all possible temperatures and chemical potentials has attracted some 



attention 



22, 23 . In 22 



one can see that, for fields greater than 10 17 G 
the critical end point and the first order coexistence region are affected. One 
observes an increase on the size of the first order segment of the transition line 
with Tcep moving to higher values as the magnetic field increases. This result 
is certainly of interest for the physics of heavy ion collisions at intermediate 
energies, where huge fields, of the order \eB\ > m\ ~ 10 18 G, can be created 
by heavy ion currents due to the spectator nucleons. The field intensity 
depends on the centrality and beam momentum so that \eB\ ~ 5m\ can be 
reached at RHIC while \eB\ w 15ml can be reached at the LHC. 
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In 24| matter is described by free gases of baryons and mesons under 



the influence of a constant magnetic field and particle rations and yields are 
calculated in such a way that the results for the temperature and chemical 
potential are chosen for the minimum \ 2 fit- The results show that, even 
for the free Fermi and Boson gas models, a strong magnetic field plays an 
important role. The inclusion of the magnetic field improves the data fit up 
to a field of the order of B = 10 19 G. For stronger magnetic fields, it becomes 
worse again. In this calculation, baryons with spin 3/2 subject to an external 
magnetic field had to be described and the degeneracy of the Landau levels 
(LL) was an important quantity. 

The Rarita-Schwinger (RS) equation Q describes spin 3/2 particles and 
has been studied widely because of its intrinsic interest and also for appli- 
cations in different research fields, and not only the above mentioned ones. 
The onus of inconsistencies represent a disadvantage to the RS equation. The 
quantization of the interacting spin 3/2 field turned out to be inconsistent 



with Lorentz covariance 191. Moreover, the wave fronts of the classical so- 



lutions of the RS equation present noncausal modes of propagation [25|, |26 



This kind of inconsistencies are old problems and several remedies have been 



suggested through the years [271]- (30 



of a nonminimal coupling [31 



including an interesting prescription 
From this last reference, it is clear that, at 
least for constant fields, the ones responsible for the noncausality problem 



25, 26 



19 , a valid, causal 



and the inconsistencies with Lorenz covariance 
Langrangian exists. 

For the reasons mentioned above, a complete solution for the RS equation 
under the influence of magnetic fields deserves a detailed revision, which we 
per form in the present work. We next follow the approach developed in 



[32J and [33| and circumvent the problems just mentioned, without having 
to rely on nonminimal couplings. To the best of our knowledge, this is the 
first time that a complete calculation of the Rarita-Schwinger equation with 
the inclusion of the magnetic interaction, including the occupation of the 
LL, is obtained in a simple way, avoiding the problems observed in older 
calculations. 



2. Formalism 



A spin-3/2 field can be described by the Rarita-Schwinger(RS) Lagrangian 
given by (with units in which H = c = 1) 

C RS {A) = r (id a T^ a V {A) - mB^(A)) y , (1) 
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where ip 11 = (if) 1 *) 7° is the adjoint vector spinor, m is the particle mass, ip v 
is the vector spinor and the matrices 1^ a v and B^ u are given by 29] : 

T, a V (A) = g, vl a + A ( 7 ^« + gfa) + B ltll a lv , (2) 

B^(A) = g^-G^ynv, (3) 

with 

B=\a* + A + \ (4) 

C = 3A 2 + 3A + 1 . (5) 

where the metric tensor is g^ u = diag (1,-1,-1,-1), 7^ are Dirac matrices 
and A is an arbitrary parameter (except that A ^ —1/2). The parameter 
A has no physical meaning and the conventional choices made in 



the literature are A = —1/3 that corresponds to the Lagrangian originally 



proposed in [9], or, A = [211 . or, for A = — 1 the Lagrangian reduces to [34 



C RS {A = -1) = -- ^ (e^ A 7 57^ p + m/ A ) ip\ (6) 

where e^ upX is the Levi-Civita symbol, 7 5 = 27°7 1 7 2 7 3 and a px = | [7 M , 7 A ]. 
The freedom represented by the parameter A reflects invariance under rota- 
tions mixing the two spin-l/2 + and l/2~ sectors residing in the RS repre- 



sentation space besides spin-3/2 [35(- 36]. By applying the Euler-Lagrange 



equations to (EJ), the equation of motion are given by 

(e^ A 757^ + m^ A )^A = (7) 

while (JZJ) can be rewritten in a simpler form plus supplementary conditions 

(iy^ _ m )^ v = , (8) 

7^ = 0, (9) 

d^ = 0. (10) 

These last results clearly show that the equations of motion (j7]) can be rewrit- 
ten as a Dirac equation for the vector spinor plus supplementary conditions. 
Therefore, the solution for this class of systems consists in searching for solu- 
tions for four Dirac-type equations. It should not be forgotten, though, that 
the equations are not indeed independent because they have to satisfy the 
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constraints OH]) and f llOp . These constraints are necessary to eliminate the re- 
dundant components of the vector spinor ?/> M . These constraints, in general, 
do not hold in the presence of interactions. Nevertheless, it is possible to 
derive the necessary number of subsidiary conditions for at least a selected 



class of interactions [37|]. Furthermore, by multiplying equation flBj) from the 
left by 7^ and using equation ([§]), we obtain equation ffTDl) . Therefore the 
"gauge" condition (IT01) is not an independent relation. The equations ([S]) and 



are known as the RS equation [38 



The RS equation can also be derived from the Lagrangian 0] 
C RS = + im) + l^d ■ + ~ ^7 M {$ -im)(f>, (11) 

where = 7^*9^, <p = T^V^ an d d-ip = d^ip^. In this case, the Euler-Lagrange 
equations of motion yield 

{$ + im) ^ - - {d^cj) + ^d-i)) + \ rf ($ - im) <j) = . (12) 

By contracting equation ( fT2"l) with <9 M and 7 M one obtains the supplementary 
conditions ([9]) and (flOl) . which allow equation ( fl2l) to reduce to the simpler 
form for the RS equation. Due to simplicity once both equations are equiv- 
alent, we shall henceforth choose the equations ([8]) and Q to describe a 
spin-3/2 particle. 

Now let us consider the system embedded in an electromagnetic field. 
This is implemented via minimal coupling prescription, which consists in the 
replacement 

d /t ^D fl = d /t + iqA /t , (13) 

where is the /i-th component of the four-vector potential that represents 
an external electromagnetic field and the charge q = e q \q\, where e q = +1(— 1) 
corresponds to a particle with positive (negative) charge, respectively. So, 
the RS equation with an electromagnetic field becomes 

{iYDf, - m)ip u = with 7^ = . (14) 

According to the previous discussion on the change of the subsidiary condi- 
tion, we can calculate the new condition by multiplying the first equation in 
(|T4"]) from the left by 7^ and using the second equation in (Tl4"|) . leading to 

D^ = 0, (15) 
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which is in agreement with the minimal coupling prescription given in equa- 
tion (JTUJ). This last result clearly shows that the minimal coupling prescrip- 
tion is mathematically consistent, in contrast to the introduction of a Lorentz 
scalar potential, which gets to a contradiction if we do^the same in the RS 

To achieve 



25. 26 



equation written in the form (joj), as addressed in 
the solutions of the RS equation embedded in an electromagnetic field, we 
follow the approach used in obtaining the solutions of the Dirac equation in 
32]. A similar approach was also used in jiif . 



Because of the reasons mentioned in the Introduction, we focus our at- 
tention on considering the external vector field to be a constant magnetic 
field. Assuming the gauge = d^B, i.e. = (0, 0, xB, 0) that yields a 
constant magnetic field transverse to the xy-plane, we have V • A = and 
B = V x A = Be^. In this manner, the covariant derivative becomes 



D, 



da - ie q \q\Bx5 i 



fi2- 



(16) 



As we have chosen the constant magnetic field to be transverse to the xy- 
plane, we can apply the method of the separation of variables in the vector 
spinor r/>^, which can be rewritten in the form 



cEt 



f^(x)v 



-ieEt+iep y y+iep z z 



(17) 



where we define E as positive, and e = +1(— 1) corresponds to the states of 
positive (negative) energy, respectively. 

Taking advantage of the solutions for the Dirac equation, we can consider 
the solution for e = +1 as 



f(e) 
J u 



( fA x ) \ 




fA x ) 













and for e n 



-1 as 



f(e) 
J u 



( fA x ) \ 












\ M x ) I 





C£?«„-l(0 



(19) 
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where 



1 



■H n (£)e 



-e/2 



(20) 



v / 7r 1 /22«n! 

is the solution for the one-dimensional harmonic oscillator, where energy 



takes the form 



E = y/p 2 z + m 2 + 2n\q\B 



and the coefficients Cji satisfy the relation 



(21) 



(e) 
/u4 



\ 



1 



ep z 

eE+rn 
eE+m 



J 





1 

eE+m 
\ eE+m 



c 



111 ■ 



(22) 



We have now to consider the subsidiary conditions. We can define the oper- 
ators 



Ox = i(ep y - e q \q\Bx + 



dx' 



i(\q\BY/ 2 (-eJ + 



2 

noticing that 
and 



-ep y + e q \q\Bx + 



dx' 



d_, 
d 



i(\q\B) 1/2 (e g Z+^), 



OlV r , 



2 v n 



-ip n+ iv n+1 

iPnVn-l 
ipnVn-1 

-ip n+1 v n+1 



if 
if 

if 
if 



"I, 

= 1, 
= -1, 



(23) 
(24) 

(25) 
(26) 



where p n = yj2n\q\B. Hence, the condition (115p (ie. D^ijj^ = 0) becomes 



eE 



+ 



Q)3 V n q (0 

(<0 



1 

+ 2 



Pn q 
Pm q 



•/ v H Cyl ) V n,_ 

-\{ie q C x2 — C y2 )v m 



-i(ie q C xl 

.(e) 



Pn q + l(itqC x3 — C 



\Pn q + l{itqC x l — C, 



( Vr, 
Prr, 
Vr, 

\Prr, 



^ £ lVn g -l(0\ 
'y3 ) V n 



(0 



+ epz 



C$)Vm q -l(0/ 



y3) V n q +l(0 
i e l)Vm q +l(0J 

C Z 2 v m q (0 
\CSvm q (0j 



(27) 



where we define n q = n (n q = n — 1) and m q 
charge is positive (negative), respectively. 

Finally, from the restriction (|9])(ie. 7 M Vv = 



/ 



"01 

r {e) v 



\ 



\ -cfiWo y 



/ 



+ 



\ 



\ 



t(e) 
x4 

y 



/ 



= n — 1 (m q ■■ 
0), we obtain 

(0\ / 



n. 



when the 



^y4 v m q 



a 



CO 



V 



°3/l y n. 



(0 
(0 

,(o y 



2:4 V " 



\ 



\ 



(0 

(28) 

From the system of equations (I27p and (1281) . we conclude that for n > 2 
only four coefficients are independent meaning that the states of positive or 
negative energy are quadruply degenerate. On the other hand, for n — 1 



(considering e q 



-1) from the system of equations fl27j) and fl28|) (t>_i 



0), we obtain that only three coefficients are independent meaning that the 
energy state for n — 1 is triply degenerate. Finally, for n = (considering 

0, and in this case 



from equation (TE?]) we obtain that Cj^ 



from (12T|) and (12"8"|) (f _ 2 = f-i = 0) we conclude that only two coefficients are 
independent, meaning that the energy state for n = is doubly degenerate. 
This degeneracy of the energy state can be associated with a spin label, as 
well as 

1 



n 



I e + 

2 q 2 



(29) 



where I = 0, 1, 2, . . . and s = ±1, ±3. Therefore, the degeneracy factor for 
spin 3/2 particles are given by 



9n — 4 — o n i — 2o n o . 



(30) 



Explicitly, we obtain g = 2 (n = 0) for Z = 0, s = —1 and I = 1, s = —3. 
We obtain g% = 3 (n = 1) for Z = 0, s — +1; / + 1, s = — 1 and Z = 2, s = —3. 
And finally, we obtain g n = 2s + 1 = 4 (n > 2) for s = 3/2. At this stage, it is 
useful to mention that independent coefficients may be completely calculated 
following an analogue procedure to the one described in the section 5.2 of 
the Ref. 4l|], i.e, using an appropriate helicity operator for spin 3/2 particles 



42j , but these calculations are beyond the scope of this work. 



0. 



3. Final discussions 

In the present paper, we revisit the Rarita-Schwinger equation and include 
the electromagnetic interaction through the minimum coupling scheme. We 
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assume a gauge where the magnetic field transverse to the xy-plane is con- 
stant. This choice leads to Dirac-type solutions, which can consistently avoid 
noncausality problems. We must observe that the energy given in eq.( l2ip may 
assume complex values for n > 1 if B > m 2 /2n\q\, but this is a very strong 
magnetic field (|el?| ~ 10 20 G), much stronger than the estimates for the 
typical values of the fields produced in heavy ion collisions 24j for example. 
Moreover, we can clearly see that the occupation of the Landau levels by 
spin 3/2 particles is quite different from the occupation of spin 1/2 particles. 
While systems composed of spin 1/2 particles subject to magnetic fields can 
accommodate only one particle in the lowest level and two particles in the 
other levels, systems with spin 3/2 particles can accommodate 2 particles in 
the first level, 3 particles in the second one and 4 particles in the remaining 
levels. This fact has a great influence in studies involving equations of state, 
density of particles, scattering amplitudes and other calculations of interest 
for nuclear astrophysics and physics of heavy ion collisions mentioned in the 
Introduction. 
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